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The theory of magnetically accelerated outflows and jets from accretion disks is 
reviewed at an introductory level, with special attention to problem areas like the 
launching conditions of the flow at the disk surface, stability of the magnetic field, and 
coUimation mechanisms. This text will appear in R.A.M.J. Wijers, M.B. Davies and 
C.A. Tout, eds., Physical processes in Binary Stars, Kluwer Dordrecht, 1996 (NATO 
ASI series). 



1 Introduction: the case for magnetic accelera- 
tion 

Narrow, high speed outflows (jets) and less weU coUimated 'bipolar' outflows are 
observed from very different cosmical objects, ranging from protostars in the 
solar neighborhood, to galactic X-ray binaries, to the nuclei of active galaxies. 
The magnetic acceleration mechanism for outflows from accretion disks has 
gained significant popularity as an explanation for each of these forms of outflow. 
The model can account for high speeds (for example, Lorentz factors of 10 or 
higher in AGN and galactic black hole binaries), high degrees of collimation, and 
large momentum fluxes. Though other processes can also, to varying degrees, 
account for these properties (e.g. Blandford 1993), the magnetic model combines 
them in a natural way. While some processes are very different in protostellar 
disks and the AGN or X-ray binary disks, in particular those that produce the 
observed radiation, the physics of the magnetic acceleration model is to a large 
degree independent of these. Progress made by development of the model for 
explaining observations in one area is therefore likely to have impact for the 
interpretation of other outflows as well. 

In spite of the rather general applicability of the magnetic mechanism, it is 
unlikely that it is involved in all cases. There are a number of bipolar looking 
objects in the sky where the explanation may well be purely hydrodynamical, 
such as the outflows in rj Carinae, and in particular the planetary nebulae (Icke 
et al. 1992). 

The protostellar outflows play a key role in supporting the magnetic/centri- 
fugal model. The momentum flux in these objects can be measured from flow 
speeds and inferred mass densities, and in many cases turns out to be much 
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larger than can be accounted for by the nearest competing mechanism, radiation 
pressure from the central star (for a review see Konigl and Ruden, 1993). The 
protostellar outflows are also the ones in which observations are most likely, 
in the near future, to reveal their inner workings by direct imaging, as shown 
in table |l|. See the contributions elsewhere in this volume for more about this 
subject. 



2 Presence of jets and outflows in binaries 

Jets are now known from all classes of X-ray binaries, i.e. mass transfering 
binaries (Lewin et al. 1995) in which the primary is a black hole or neutron 
star. (Ignoring here the 'supersoft' X-ray binaries in which the primaries appear 
to be white dwarfs). Among the massive X-ray binaries (HMXB) in which the 
companion is an early type star, there are the well known jets of SS 433, and the 
radio jets from Cyg X-3 (Strom, van Paradijs and van der Klis, 1989) and the 
galactic center source 1E140.7 - 2942 (Mirabel et al. 1992). Among the low mass 
X-ray binaries (LMXB) with neutron star primaries, a jet is known only from 
Cir X-1 (Stewart et al. 1993). Until recently no jets were known from LMXB 
with black hole candidate primaries, but this has changed with the discovery 
of superluminal jets in the variable source GRS 1915-f 105 (Mirabel 1994), and 
the transient GRO 1655-40 (= X-ray Nova Sco 1994), (Hjellming and Rupen 
1995). The physics in the inner disks of these objects must be rather similar to 
that in the central engines of AGN (e.g. Begelman, Blandford and Rees 1984), 
and it is pleasing to see that they can produce very similar jets, though on a 
much smaller scale (for more observational similarities between galactic black 
hole candidates and AGN see Sunyaev et al. 1991 and references therein). 

Outflows without evidence for jets exist in Cataclysmic Variables (CV: bi- 
naries transfering mass from a main sequence star to a white dwarf, see Hack 
and La Dous, 1993). P-Cygni proflles indicating mass loss are seen in UV lines 
in Dwarf Novae (DN) when in outburst (e.g. Drew and Verbunt, 1988). No 
evidence for outflow is known for DN in quiescence. In UX Uma systems (CV 
with steady mass transfer), evidence for mass loss comes from single-peaked line 
profiles (in contrast with the classical double-peaked profiles of accretion disks). 



Table 1: Angular size of the accelerating region (assumed to be 100 times the 
typical size rg of the inner disk) for jet-producing objects. 

inner disk distance angular scale (") 
ro D lOOro/D 

nearby protostar 3Rq 500 pc 0.003 

nearby active galactic nucleus 100 AU 10 Mpc 0.001 
galactic black hole candidate 100 km 2 kpc 3 10~* 
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and the presence of additional 'uneclipsed light' in eclipsing systems ('SW Sex 
syndrome', Thorstensen et al. 1991). It is fair to say that no collimated outflow 
has yet been observed from a CV. One might wonder if this has something to do 
with the fact that the primaries are white dwarfs, but the case of R Aqr shows 
that this cannot be the case. R Aqr is a white dwarf with a Mira type giant 
companion in a 44 yr orbit; mass transfer occurs because of the dense stellar 
wind from the Mira. It has a jet which is visible at optical as well as radio 
wavelengths (Burgarella and Paresce 1992, Dougherty et al. 1995). 

It is still somewhat puzzling that among the bright LMXB containing neu- 
tron stars there is only one case with a jet (Cir X-1). Also, it is not clear why 
jets have been found only in a few of the relatively frequent (1-2 yr^^) bright 
X-ray transients. These systems are believed to be all rather similar. This indi- 
cates that there may be additional factors determining the production of jets, 
factors for which no observational counterparts have been identified so far. 

Protostellar binaries are a bit of a different class of objects in this context, 
since their disks are not fed by mass transfer. The effects of a binary companion 
on a protostellar disk, however, may well be relevant for the ability of the object 
to produce a jet. I return to this question in the last section. 

3 Physics of magnetic acceleration: heuristics 

Magnetized winds and jets can be produced by rotating objects which, for one 
reason or the other, have a magnetic field anchored in them. The importance 
of such a magnetic field for the spindown of stars was realized by Schatzman 
(1962) soon after the discovery of the solar wind. Quantitative models for 
magnetized stellar winds were then developed by Weber and Davis (1967) and 
Mestel (1968). The point of view in this work was, mostly, the spindown of 
stars, but Michel (1969, 1973, see also Goldreich and Julian 1970) realized the 
importance of the mechanism for producing high speed outflows, and formulated 
relativistic models for magnetic winds from pulsars. That strong outflows could 
also be driven magnetically by accretion disks was proposed by Bisnovatyi- 
Kogan and Ruzmaikin (1976), Blandford (1976) and Lovelace (1976). Workable 
quantitative models for such flows were first produced by Blandford and Payne 
(1982), while full numerical solutions for the steady (nonrelativistic) problem 
were first obtained by Sakurai (1985, 1987). 

In this section I review, at a heuristic level, the basic ideas and processes 
involved. The theory is introduced a bit more formally in the next section. For 
definiteness a wind generated by an accretion disk is considered here. Most of 
the discussion however, applies equally well to rotating stars. 

First, divide the disk and its surrounding space according to the relative 
importance of the magnetic field energy density, see figure |l|. The magnetic 
field to be used for producing the wind is anchored in the disk, so its energy 
density there must be less than the rotational kinetic energy in the disk. It 
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Figure 1 : Regions of force-free and non-force free magnetic field in a disk-driven 
wind. 



can exceed the thermal energy density in principle, but in any case its strength 
and configuration is determined by other forces, which provide the anchoring 
of the field. Just how large the field strength is, for any kind of observed disk, 
is still unknown since no field strengths have been measured yet^. Theoretical 
arguments allow for field strengths of the order of equipartition with the gas 
pressure, in the case of dynamo-generated fields (e.g. Hawley et al. 1995, Bran- 
denburg et al. 1995), or even larger field strengths for magnetic flux dragged in 
with the accretion flow (Spruit, Stehlc and Papaloizou 1995, Lubow and Spruit 
1995). Without being too specific, I assume that the vertical field strength, at 
the disk surface, is reasonably large, since the magnetic acceleration mechanism 
requires a field of some strength (a more specific criterion is given in section ^ . 
Outside the disk, the gas density is typically so low (assuming a cool disk) that 
the magnetic energy density is large compared with the thermal and rotational 
energies. The field in this region must therefore be force free, much like the 
solar corona (e.g. Foukal 1990). In the absence of torques acting on it, it must 
even be a potential field. As figure Q suggests, we are assuming that the field 
is of uniform polarity, over the region of the disk where the wind is generated. 
Loops of field connecting different parts of the disk cannot be excluded a priori. 
Such loops are sheared rapidly by differential rotation, giving rise to a rich and 
not very well understood complex of phenomena which is outside the scope of 
this discussion. It suffices that a certain minimal fraction of the disk's magnetic 

^An exception is the field strength (about 1 G) in the protosolar nebula inferred from 
meteorites (e.g. Cisowski and Hood 1991) 
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Figure 2: Bead-on- a- wire analogy for centrifugal acceleration by a magnetic 
field. 



flux does not loop back to the disk surface, but is open to infinity. 

As the flow is accelerated, the field strength drops due to the increasing 
distance from the disk. The acceleration effectively stops when the flow speed 
reaches the local Alfven speed in the flow. The place where this happens is 
called the Alfven surface. Thus, outside the region where the magnetic field 
dominates there is again a region where the field is not force free; in this case 
because of inertial forces (figure |l|) . 

The acceleration process is illustrated in figure Assume that the disk is 
cool, so that its rotation is close to Keplerian, and its thickness can be neglected. 
Also assume that the gas is sufficiently ionized everywhere that ideal MHD 
can be used, i.e. that gas is tied to the field fines. These assumptions are 
not essential and can be relaxed in numerical models, such as those of Konigl 
(1989). Assuming we are close to the disk surface, and the field strength large, 
the atmosphere of the disk is forced to corotate with the field lines sticking out 
of the surface. Since the Lorentz force only has components perpendicular to 
the field, the gas is free to move along the field line, under the influence of the 
other forces, like a 'bead on a wire'. At the foot point of the fleld line, the inward 
force of gravity just balances the centrifugal force, because of our assumption 
of Keplerian rotation in the disk. Along the fleld line, the centrifugal force 
increases with distance from the axis. When the component of the centrifugal 
force along the fleld line exceeds that of gravity, the gas tied to the fleld line is 
accelerated outward. 

This centrifugal process stops when the flow speed becomes comparable to 
the Alfven speed; at that point, the fleld is no longer strong enough to enforce 
corotation. 

Depending on one's preferences, the acceleration can also be described in 
purely magnetic terminology (e.g. Lovelace, Wang and Sulkanen 1987), without 
appealing to a centrifugal force. This is discussed in section p|. 
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Figure 3: Development of the azimuthal field. With each rotation of the field 
line a loop of field is added to the flow at the Alfven surface. 

Beyond the Alfven surface, the inertia of the gas causes it to lag behind 
the rotation of the field line, so that the field gets wound up. The simplest 
way to visualize this is by ignoring the rotation of the gas altogether. Then for 
each rotation of the foot point of the field line, one loop of field is added at 
the Alfven surface. As the flow carries these loops away, a spiral shaped field 
formed (figure with pitch v/Q., where v is the flow speed and Vl the rotation 
rate of the foot point. In a (nonrotating) frame comoving with the flow, one 
sees a nearly azimuthal magnetic field, its strength decreasing with time. The 
curvature force in the azimuthal field is directed towards the axis, causing the 
flow to 'coUimate', i.e. to become parallel to the rotation axis. 

A numerical example of a magnetically accelerated flow is shown in figure ^. 
Note that the rotation velocity peaks near the Alfven radius r\, which is 
at 100 times the starting distance rg in this example. At large distance, the 
rotation drops roughly in accordance with angular momentum conservation. 
The asymptotic radial velocity is larger than the rotation speed at ta by a 
factor of order unity. 

One of the attractions of the magnetic wind model is that it not only pro- 
duces outflow, but in principle can also take out the angular momentum from the 
disk, allowing it to accrete (Blandford 1976, Bisnovatyi-Kogan and Ruzmaikin 
1976). To see how effective this can be let us estimate the angular momentum 
flux carried by the wind. Since the flow corotates roughly up to the Alfven 
point, the specific angular momentum carried is of the order ^r\, hence 

Jw = I^U^rl, (1) 

where is the mass flux on a field line with foot point tq. It turns out that 
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Figure 4: Example of a magnetic wind model. Top 2 panels: rotation velocity 
(measured in an inertial frame) and radial velocity in units of the rotation speed 
at the base of the flow, as functions of distance. Alfven radius is at lOOro- Lower 
panel: field angle. (Cold WD model for /j. = 10~^, see text). 
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Figure 5: CoUimation of the flow by the curvature force of the azimuthal 
magnetic field. 

this estimate is actuahy exact (section ^) . The angular momentum that has to 
be extracted (locally at tq) from a Keplerian disk in order for it to accrete at a 
rate Ma is Ja — ^firpMa, hence 



Since this relation is exact, and ta always larger than tq, it follows that only a 
fraction of the mass flux in the disk can flow out in the wind. It is possible in 
principle, however, that the wind carries away all angular momentum that the 
disk has to loose in order to accrete, i.e. without angular momentum transport 
by viscous torques in the disk. Such a disk would, in the absence of viscous 
dissipation, be silent. This is verified by looking at the energy balance. The 
energy flux in the wind is given by flJ^, (the work done against the wind 
torque). If all the angular momentum is carried with the wind we have = E^ 
(using equation ||), where £^a = l/2(riro)^Ma is the rate of gravitational energy 
release in a Keplerian disk. Thus, if the wind carries away all the angular 
momentum, it also carries away all the accretion energy. 

3.1 Structure of a disk driven wind 

The magnetically driven wind is conveniently broken down into three conceptual 
stages or subprocesses. Near the disk surface, the wind is 'launched': details 
of disk structure and magnetic field configuration near this surface determine 
how much mass is launched into a fiow. After this, the flow can be regarded 
as ballistic, its acceleration being governed almost entirely by gravitational and 
magnetic/centrifugal forces. After the acceleration phase, which ends roughly 
at Alfven surface, the coUimation phase starts, in which the flow is deflected 
towards the axis by 'hoop stress' (figure |^). 

Since the outer radius of a disk is typically much larger than the inner 
radius, conditions can vary dramatically with distance in the disk. The wind 
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problem is, therefore, a function of distance. In the inner regions where the field 
strength and the rotation speed are large, conditions are favorable for producing 
high speed coUimatcd winds. In the outer regions, one would expect lower 
wind speeds, and probably less coUimation. The fluxes of mass and angular 
momentum from the outer regions, on the other hand, could be large compared 
to flows from the inner regions. Flows from these regions may well coexist. 
Since the magnetic field inside the Alfven radius is force free, its configuration 
is determined by a global force balance. The wind properties from adjacent 
regions in the disk are therefore coupled somewhat, through their dependence 
on the shape of the poloidal field. This aspect of the disk wind problem has not 
yet received much attention. 

In the next sections this picture is elaborated a bit more formally. In the 
process, some problem areas of current interest will be noted, relating, in par- 
ticular, to the launching and coUimation phases. 

4 Steady axisymmetric magnetic flows 

The theory of magnetically driven flows from rotating objects has been given 

in many texts (e.g. Weber and Davis 1967, Mestel 1968, Hcincmann and Olbcrt 
1978, Okamoto 1975, Blandford and Payne 1982). I repeat the basic derivation 
here, under the assumption of ideal magnetohydrodynamics, in the nonrela- 
tivistic limit. The extension to relativistic MHD can be found in Michel (1969), 
Goldreich and Juhan (1970), Bekenstein and Oron (1978), Camenzind (1987), 
for a general treatment of relativistic MHD see Lichnerowicz (1967). 

The basic assumptions made arc the MHD approximation, and that the 
flow is stationary and axisymmetry. The MHD approximation is justified if the 
rotating object actually produces an outflow of any observational significance; 
it holds if the density of charge carriers is large compared with the so-called 
Goldreich- Juhan density, Nqj = flB/{4Trce) ~ lQ~'^i}B cm~^. This hmit is of 
importance in the case of pulsar magnetospheres (Goldreich and Julian 1969), 
but is so low that it is not likely to become relevant for most of the observable 
winds and jets. Though deviations from stationarity are implied by the pro- 
duction of outward traveling 'knots' in most jets (for an example see Mirabel 
and Rodriguez 1994), the time scale for acceleration is still likely to be short 
compared to the time scale of these variations, so that stationarity is a good 
assumption during the acceleration phase. Much of the jet phenomenology is 
consistent with axisymmetry. Theoretically, nonaxisymmetric instabilities are 
likely to become important in parts of the jet; this is discussed further below. 
Keeping this in mind, we proceed with the axisymmetric case. The final assump- 
tion made is that of infinite conductivity. This can easily be relaxed (Konigl 
1989), but deviations from this approximation are important only in the case 
of protostellar jets, where the flows are so cool that the conductivity needs to 
be considered in detail. 
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Figure 6: Axisymmetric poloidal magnetic anchored in a rotating object. Mag- 
netic surfaces are labeled by the flux function ij). Rotation rate can be a function 
of ij). 

The equations for stationary ideal MHD are (e.g. Roberts, 1967) 



V X (v X B) = 0, (3) 

pv Vv = -Vp- pV$ + -^(V X B) X B, (4) 

V • (pv) = 0, (5) 

V • B = 0. (6) 



Here B, v, p, p are the magnetic field vector, the velocity, the gravitational po- 
tential, the gas pressure and the density. Using cylindrical coordinates (lu, (p, z), 
axisymmetric vectors like B and v can be written in terms of their poloidal (p) 
and toroidal (0) components: 

B = Bp + 5060, v = Vp + i;0e0, (7) 

where the poloidal components lie in the meridional {w, z) plane. 

With the assumed axisymmetry, these equations have a great deal of sym- 
metry, so that they can, in part, be reduced to algebraic eriuations. This is done 
as follows. Due to axysimmetry and divB = 0, Bp can be written in terms of a 
flux function tjj: 

Bp = — X e^, (8) 

zu 

i.e. Bz = 1/w dip/dw, Bj. = —1/w dip/dz. The flux function (also called 
stream function or vector potential) is constant along field lines: 

B • W = 0. (9) 
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Thus, ip can be read as a label numbering the field Hues of the poloidal field 
(figure ^). It plays an important role in this numbering capacity, since there 
will be several scalar fields a with the property Bp • Va = 0. Such fields have 
their gradients parallel to that of V', and therefore are functions of ip only, and 
are also constant on field lines. From (equation ^ we have B x v = V/, 
where / is an axisymmetric scalar. Writing this out into poloidal and toroidal 
components: 

Vp X Bp + v^e^ X Bp + B^Vp x = V/. (10) 
The toroidal component of this is 



'p 



X Bp = 0, or Vp = K(n7, z)Bp. (11) 



Thus, the poloidal velocity is parallel to the poloidal magnetic field, a conse- 
quence of the infinite conductivity assumed. With (0), the dot product of (|o|) 
with Bp yields Bp • V/ = 0, so that / is not just a scalar, but a function of the 
field line number ijj only. With this, (|l^) yields 

KB^^ruf'i^j), (12) 

where /' ~ df/dijj. With ([III), (||), the continuity equation (||) yields Bp • 
V(/9k) = 0, so that 

pVp/Bp = pn^ ri(jp), (13) 

for some function 77. This equation simply states that the mass flux density, per 
unit of poloidal magnetic flux, is constant along a field line: each field line has 
its own mass flux, in this sense. 

Now find the point along a poloidal field line inside the rotating object, 
on which — 0, and call this the 'foot point' of the field line. If the object, 
together with Bp, is symmetric about the equator, this point is on the equatorial 
planeF. Let J7 = v^q/wq be the rotation rate of this foot point. Then we find 
from (12|) that f'{ip) = O. Loosely, we can call this the 'rotation rate of the field 
line'. Remember, however, that the plasma rotation rate, f2p — v^/w, is not 
constant along a field line. It cannot be, because the plasma corotates with the 
foot point only as long as the field is strong enough to dominate over the plasma; 
the rotation starts lagging behind at larger distances from the object where the 
field is weaker. Let v' be the fiow velocity measured in a frame rotating with 
the angular frequency Q.: 

v' = V - wVt{il;)e^ (14) 
Then (|2|) and ( |ll| ) can be combined into 

v' = kB. (15) 

^Such symmetry simplifies the visualization, but is not necessary. All poloidal field lines 
are 'bent back' at the object's surface by the wind torque in the same direction (lagging 
behind the rotation), so that is of opposite signs at two successive crossings of the field 
line through the object's surface. Hence on all poloidal field lines crossing the surface the 
toroidal component changes sign at some point inside the object. 
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In other words, in a frame corotating with the foot point of a field hne, the flow 
is everywhere parallel to the magnetic field (this frame can be different for each 
poloidal field line, if the object rotates differentially). Whereas the poloidal flow 
component is parallel to the poloidal field component in both the rotating and 
the stationary frames, the total velocity is parallel to the total B-field only in a 
corotating frame. 

Next, consider the equation of motion (in a non-rotating frame). Using 
axisymmetry and the identity (V x B) x B = — + (B • V)B, the toroidal 
component is 

p(v-Vv)^-i-(B.VB)^. (16) 

With the vector relation (a • Vb)^ — a • 'S/{-wb^) /w, (^ij) and (p^), this can be 
written as 

Bp • ^{pnvjv^) = —Bp • V(roB^). (17) 
47r 

This can be integrated, yielding 

-^[pvpwv^ - ^B^Bp) = i]L, (18) 

where L is a a function of ■0 only. Hence 

^{v^ - -^B^) = L{iP). (19) 

The first term in (^) is the flux of angular momentum by the poloidal flow, the 
second the magnetic torque. Thus, ( [iq ) expresses that the total flux of angular 
momentum, per unit of poloidal magnetic flux, is constant along each field line. 
The total angular momentum flux measured per unit of mass flowing along the 
field line is L. 

Before proceeding with the equation of motion, we use (^|) to eliminate B^ 
from (|l|). This yields 

L - nw^ 

"^'""^^ ..[i-i/(Wp)] - 

The denominator on the RHS vanishes when AnpVp/Bp = 1, or Wp = UAp, where 
i^Ap = Bp/ (47r/9)^/^ is the poloidal Alfvcn speed. The location along the field line 
where this happens is called the Alfven point, because an axisymmetric Alfven 
wave propagates along a magnetic surface at the speed WAp (independent of the 
value of B^). If the flow is to be accelerated to values beyond the local poloidal 



Alfven speed, the numerator in (20) also has to vanish at the Alfven point. If 
we denote the physical variables at the Alfven point with a subscript a, this 
condition yields 

L = (21) 
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The interpretation of the singularity at wk is similar to that of the sonic point 
in the transition from subsonic to supersonic flow in ordinary hydrodynamics 
(jet nozzle, stellar winds). Eq. ( pl| ) has an important physical interpretation. 

Eq. (^l|) has an important physical interpretation. If the flow were to 
corotate with the field up to vox (which actually it does not, see above) it 
would, at that point, have the specific angular momentum given by (^l|). Since 
L is the specific angular momentum flux in the wind (including the magnetic 
torque!), it is as if the flow were kept rigidly corotating out to toa, and then 
released without further magnetic torques. As far as the angular momentum 
flux is concerned, the flow 'effectively corotates' out to wa- 

As in ordinary dissipationless hydrodynamics, the equation of motion in the 
direction of the flow can be integrated in terms of a Bernoulli function, provided 
that the energy equation is sufficiently 'simple'. In practice this is the case if 
polytropic or isothermal equations of state can be used as approximations: 

P = Kp^, (22) 

where /C, 7 are constants. The isothermal case is obtained with 7 = 1. These 
are not necessarily very good approximations to the real situation, since energy 
dissipation and cooling processes usually are present. These play an important 



role in the early stages of the acceleration, see section 5.3. Once the flow has 
been accelerated beyond the speed of sound, however, the dynamics does not 
depend much on the temperature of the gas any more. 

In a corotating frame (corotating with the foot point of the field line, cf. the 
discussion above), the equation of motion becomes: 

pV ■ Vv' = -Vp + — (V X B) X B - pV$ + pn^w + 2pV x n, (23) 
47r 

where w = we^ and $ is the gravitational potential: 

GM 

$ = . (24) 

r 

Let s be a unit vector parallel to the magnetic field. Taking the dot product 
of (|23| ) with s, we get the component of the equation of motion parallel to the 
field lines. I denote the derivative along the field line by 

ds=s-V. (25) 

Then the centrifugal term fi^s • m can be written as dsi^'^vu'^ /2), since J7 is a 
constant along a field line (remember it is the foot point rotation rate, not the 
local fluid rotation) . The Coriolis term disappears because it is perpendicular to 
V, which is parallel to B in the rotating frame. The magnetic term disappears 
because the Lorentz force is perpendicular to B. Thus (|23| ) yields 

d,^v'^ ^~-^dsP~ds^ + ds{^n^w^). (26) 
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With the polytropic equation of state, the thermal term can be written as 

ia,p = 9,(fcc^), (27) 



where fc is a factor of order unity: 



(7 ^ 1), (28) 



7-1 

= ln(p) (7-1), 

and 

Cs = {vlpfl^ (29) 

is the isothermal sound speed. The equation can now be integrated; with the 
definition t;', this yields: 

+ ^(«0 - ^^f + H + * - \^^^'' - (30) 



alternatively: 



2 P 2' ^ ' 2 

1 



^-w -v^l^ro + Zcc, +$ = £;(V'). (31) 

The integration constant E depends on the field line label -0 only. This is called 
the Bernoulli equation. It states that, in the rotating frame, the sum of kinetic, 
thermal, gravitational and a 'centrifugal energy' is constant along a field line. A 
look at the terms in this equation shows the basics of centrifugal acceleration. 
Assume that the field is strong enough to enforce approximate corotation, so 
that « flm. Assume that we are looking at a flow which has already been 
accelerated to supersonic speeds, Vp ^ c, so that the thermal term can be 
ignored compared with the first term. The kinetic energy then increases 
with distance zu from the axis due to the rapid decrease of the centrifugal term. 
This is offset by the increase of the gravitational potential but eventually the 
centrifugal term dominates because $ reaches a constant value at infinity. More 
precisely, the condition for outward acceleration, at any point on the field line, 
and still ignoring thermal effects, is: 

9,($-ir!W)<0 (T = 0), (32) 

where the arc length s is taken to increase in the direction of va. The thermal 
term adds to the acceleration: as the flow expands outward, thermal energy 
is converted into kinetic energy, in the same way as in a jet nozzle. In slowly 
rotating stars, this is the dominant process driving the stellar wind, and one 
has a thermally driven wind (Parker 1963, for an introduction see Foukal 1990). 
In more rapidly rotating stars, and in our case of jets produced by disks, the 
thermal energy plays a role only in the initial stages, and most of the actual 
acceleration is centrifugal. The conditions for launching a wind from a disk are 



discussed further in section 5.2 
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4.1 Acceleration: centrifugal or magnetic? 

This discussion above suggests that the centrifugal force could accelerate the 
flow indefinitely, but this is an artefact of our assumption of corotation. When 
the field becomes weak with distance, the azimuthal velocity starts lagging be- 
hind Oci7, and then becomes small compared with Qw. The second term in 
equation |3^ then nearly cancels the centrifugal term. How much acceleration 
still remains depends on the details of how fast decreases, and the other 
equations have to be used as well to determine this. 

By working in the rotating frame there is no contribution from magnetic 
forces. This may seem strange, since it is ultimately the magnetic forces that 
transmit the rotational energy of the object in which they are anchored to 
the flow. That the acceleration can also be regarded as magnetic is seen by 
considering the equation of motion in an inertial frame. We want to know 
how the poloidal velocity is accelerated by the magnetic field. The poloidal 
component of the Lorentz force is 

Fp = T-[(V X B) X B]p = -5-[(V X B)p x B^ + (V x B)^ x Bp], (33) 

where p and ^ are the poloidal and toroidal components of vectors, as defined 
above. The second term is perpendicular to Bp and Vp, so does not contribute 
to acceleration. Thus the accelerating force is the first term on the right in (^3|), 
which can be written as 

-^(VxB0)xB^. (34) 

The condition that the poloidal flow is accelerated is then Vp • (V x B^) x B^ > 0. 
With (0) this can be written as: 

This shows that it is the pressure gradient (first term) and tension force (sec- 
ond term) of the toroidal field which determine the acceleration. For a net 
outward acceleration to occur, has to decrease outward along the field line 
sufficiently rapidly to overcome the tension force, which is directed towards 
the axis. Whether this is actually the case can not be determined from this 
argument, since one has to solve the full problem to find B^. 

We have derived both acceleration conditions ( ^2|) and ( ^5|) from the same 
equation of motion, hence the magnetic and centrifugal points of view are equiva- 
lent. This can be verified by deriving the Bernoulli equation in an inertial frame. 
The component of the equation of motion parallel to v then has a magnetic term 
V • (V X B) X B/(47r) instead of the centrifugal term. Using (|3|), (|l^), (|l8|), 
the magnetic fields in this term can be replaced by velocities. The end result is 



equation 31 
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The situation at hand determines which of these views is more appropriate. 
In regions where the field is strong enough to enforce approximate corotation, 
the centrifugal view is appropriate. When corotation is not a good approxi- 
mation, the acceleration is more conveniently viewed as due to the magnetic 
pressure of the azimuthal field. Corotation is usually a good approximation up 
to the Alfven radius (with a significant exception, see section ^ below). Beyond 
the Alfven radius, the field lines stop corotating, and instead are rapidly wound 
up into a nearly toroidal field. Here, some residual acceleration by the gradient 
of takes place. 



5 Acceleration in a fixed poloidal field 

To the extent that the field above the disk can be approximated by a potential 
field, it depends only on the distribution of its 'sources' on the disk, namely the 
normal field component i3z(n7, z ~ 0). The accelerating flow and the toroidal 
field which develops in it, however, exert forces which distort the poloidal field. 
A full solution of the problem therefore requires solving the equation of motion in 
the direction perperdicular to the poloidal field fines, the 'cross-field' equation. 
This is considered below (section A convenient approximation for the wind 
problem, however, is to consider the poloidal field as fixed and given. The 
toroidal field (which is responsible for the acceleration), is left free, to be solved 
for. This approximation is good for dealing with the launching and acceleration 
aspects of the problem in the case when the Alfven radius is at a large distance, 
since most of the acceleration then takes place in the magnetically dominated 
region. It obviously breaks down whereever the toroidal field dominates over 
the poloidal component. Thus coUimation of the flow by the 'hoop stress' of 
cannot be dealt with in this approximation. It also fails, in the entire domain, 
in the high-mass loss regime discussed in section |^. 

Having dispensed with the cross-field equation by the fixed-poloidal-field 
approximation, the flow can be solved on each poloidal field line separately. 
The solution of this problem is determined by eqs. (|30|), (20), ( |l3|) above. On 



each of the field lines, ry, wa and E are constants, still to be found from the 
solution. It turns out (Sakurai, 1985) that the problem is visualized conveniently 
in a space in which the coordinates are s, the arc length along a poloidal field 
line, and p, the gas density. The Bernoulli equation can then be read as an 
algebraic equation specifying a relation beteen p and s. This relation is the 
solution of the problem. To see this, note that for each field line, Bp is a known 
function of s, and tn is a known function of s through the known shape of the 
poloidal field lines. Hence with (|l^) the first term in ( ^ ) is of the form fp/2 = 
/i(s, p; rj). In this notation, the semicolon separates the the coordinates s, p from 
the parameters ri,E,WA- With the rotation rate of each field line specified, 
eqs. (p0|),(|2l|) show that the second term is of the form v'^ /2 — /2(s, p; -cua). 
The gravitational and centrifugal terms are functions of s only, and the thermal 
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term is of the form /3(p; K) in the polytropic case, or /3(p; Cs) in the isothermal 
case. The conditions at the surface of the disk have to be known to solve the 
problem, so we may assume that the values po,Po of the pressure and density at 
some point near the base (s = 0) of the flow are known. One of these 2 values 
then determines K or Cg directly while the other, po say, fixes (through the 
solution of the problem) a relation between the unknown constants (r],WA,E). 
Thus, for a solution of the problem, two more conditions are needed to specify 
all three constants. These are two critical point conditions, which appear as 
follows. 

5.1 Critical points 

Writing the Bernoulli equation in the form 



the solution curve p{s) can be regarded as a contour line of the Bernoulli function 
H in the s, p plane. The astrophysically relevant solutions start at a very high 
density near the disk surface, and decrease to vanishing density at infinity. We 
are therefore interested in unbroken contours of H which cover the entire p 
coordinate. H has 'mountain ranges' however, and brief inspection will show 
that these mountain ranges can be crossed by a level countour only through 
mountain passes. These mountain passes are critical points of the saddle type. 
As we shall see, there are two of these points, and the solution has to cross 
both. The elevation of H at one of the points determines the value of E. 
For the solution also to cross the other point, this point must have the same 
elevation. This will in general be the case only for certain combinations of the 
parameters ry, wa- Together with the given value of po, the two critical points 
thus determine the unknown parameters 77, wa and E and the problem is solved. 

Notice that there is no need to solve any differential equations, the whole 
magnetic wind problem is algebraic, for any configuration of the poloidal field 
(as long as it is prescribed in advance). 

It remains to be shown that there are two relevant critical points. At a 
critical point we have 



where ^Ap, I'A^ are the Alfven speeds based on the poloidal and toroidal field 
strengths, respectively: 



H{s,p;r],WA) = E, 



(36) 



dH/ds 0, dH/dp = 0. 



(37) 




(38) 



(39) 
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To interpret the expression for dH/dp, note that the dispersion relation for 
magnetosonic waves in a homogeneous medium is (e.g. Plumpton and Ferraro, 
1966): 

- {cl +v\)+ clv\ cos^ 6* = 0, (40) 

where 

u = uj/k (41) 

is the wave speed|^ and 9 the angle between the magnetic field vector and the 
direction of the wave vector k. Thus, dH/dp = when the poloidal velocity 
equals the speed of a magnetosonic wave propagating parallel to the poloidal 
flow (so that cos 6* = B^/B). Thus ( p^ ) can be written as: 

dp vi - <p 

Where Wsp and wtp are the solutions of ( ^0| ) for cos 9 — Bp/B. Critical points 
therefore occur when the flow just balances the speed of a magnetosonic wave 
propagating opposite to the flow. They are called the slow mode critical point 
and fast mode critical point, or slow and fast point, for short. In addition to 
these critical points, there is a singular point of a different kind at |tip| = fAp- 
Note that only the slow and fast mode critical points yield constraining rela- 
tions for the solution, the Alfvcn point does not yield an additional constraint. 
Through its appearance in the denominator in (|4|), the Alfven point is a node 
rather than a saddle point; all solutions which pass through the slow and fast 
points also pass through the Alfven point. This is because a critical point con- 
dition has already been applied at the Alfven point in deriving (|2l|) . The Alfven 
point, however, plays a new role as a critical point when the cross-field balance 
is considered (section H) . 

The practical problem of determining the location of the critical points and 
computing the full solution depends a bit on the character of the poloidal field 
specified. A simple case is the Weber and Davis (1967) model. This is discussed 
further in section below. 



5.2 Multiple critical points 

Since the position of the critical points depends on the geometry of the poloidal 
field configuration, one may wonder if there could not be more than just 2 
critical points. H could have an additional mountain range such that there are 
two slow points and a fast point, for example. Since each critical point adds a 
condition that has to be satisfied by the fiow, and the problem is just closed with 
two critical point conditions, one might think such multiplicity is excluded. The 
fiow, however, has a way of generating additional degrees of freedom. Except in 
carefully construed cases, a supersonic flow develops a shock at a location where 

•^magnetosonic waves, though anisotropic, arc nondispcrsive. 
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it is forced to decelerate. Thus, if the shape of the pofoidal field is sufficiently 
complex that it forces the flow to decelerate somewhere, a shock is formed near 
that location. This can happen if the flow diverges sufficiently rapidly due 
to a decrease in poloidal field strength, or if the path of the field line is such 
that it takes the flow up and down theeffective potential (gravitational plus 
centrifugal) more than once. The Hugoniot conditions will flx the properties of 
this shock, but the position of the shock is determined only if one more condition 
is imposed. The flow can pass through several critical points in such a way that 
each additional critical point is associated with a shock. After each shock, the 
flow is reaccelerated and passes though a new critical point. The regularity 
condition at the additional critical point determines the position of the shock. 
There are two kinds of shocks in magnetohydrodynamics, named slow and fast 
shocks (e.g. Jeffrey and Taniuti, 1964), so that the formation of shocks can in 
principle take place in association with additional slow critical points, as well as 
with additional fast points. An example of an additional slow point in a realistic 
disk field geometry is given in Cao and Spruit (1994). Cases with multiple fast 
points are discussed in Heyvaerts and Norman (1989). 

5.3 Launching of the wind: the sonic point 

If the flow is to be strongly accelerated, the Alfvcn point must be at a large 
distance, and this requires the field to be strong enough to enforce corotation 
out to that large distance. In this case, it is a good approximation to assume 
that the Alfven speed is large compared to the sound speed, at the base of the 
acceleration region. The slow mode speed, measured at the slow mode critical 
point, is then close to the sound speed, and the slow mode has the character of 
a sound wave guided along the fleld line. For this reason, the slow point is also 
called the sonic point. Without loss of physical generality, I equate the slow 
mode speed to the sound speed for the rest of this section. 

The importance of the sonic point for the wind problem is that it regulates 
the mass flux on the field line; it governs how much mass is 'launched' into the 
accelerating region. At the sonic point, Vp = Cg, and the mass flux (per unit 
area) is 

rh = (pcs)c, (43) 

where the index c means evaluation at the sonic point. The sound speed can be 
assumed to be known, either explicitly if an isothermal model is used, in general 
by the energy balance model used. The mass flux is then known if the density at 
the sonic point is known. At the sonic point, the pressure balance along the flow 
is afii'ected by hydrodynamic forces, but in the subsonic region before the sonic 
point their influence is small. As a fair (order of magnitude) approximation, 
we can take the pressure distribution to be hydrostatic between the foot point 
(index o) and the sonic point, and supersonic beyond. An estimate of the mass 
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Figure 7: Variation of the effective potential with arc length along field lines 
of different inclinations. 



flux is then obtained as 

TO « poCs exp[-($ec - *eo)/Csc] (44) 

where the temperature has been approximated as constant, and $c is the effec- 
tive potential including the centrifugal term, 

<fc(.) = -^- (45) 
r[s) 2 

To complete the estimate, we need a value for the potential at the sonic point. 
Assume that the potential has a maximum, measured along the field line; first 
increasing due to the gravity of the central object, and then decreasing due to 
the centrifugal force. Along the increasing part, the density is stratified nearly 
hydrostatically. Approaching the maximum, the mass starts flowing when the 
thermal energy in the gas becomes comparable to the distance to the top of the 
potential barrier. At low temperatures, this happens close to the maximum of 
the potential. An approximate mass flux is therefore found from (^) by taking 
for <i>ec the maximum of (f>o- For higher temperatures, the sonic point occurs 
somewhat before the maximum of the potential. In the absence of rotation the 
potential is due to gravity alone, and its maximum is at infinity. In this case, 
the sonic point occurs roughly at the point where the thermal energy is equal to 
the depth of the gravitational potential. Such a flow is a thermally driven wind, 
like the Sun's (e.g. Foukal, 1990). At low temperatures, the mass flux is very 
sensitive to the height of the potential barrier, since it comes in exponentially. 

For an understanding of the launching of the wind from a disk, we have 
to look more closely at the variation of the effective potential near the base of 
the flow. In flgure ^ the variation of the potential is sketched for three paths 
starting at the same point at the midplane of the disk. If the path is vertical (B), 



20 



the potential increases monotonically, there is no maximum sufficiently close to 
the disk surface, and at best a feeble thermally driven wind is possible. At an 
intermediate inclination (C), there is a maximum near the disk surface, and 
conditions for launching a wind can be good, depending on the temperature of 
the gas. Along a path close to the surface (A), the effective potential decreases 
monotonically. In this case, the wind can start right from the disk surface. The 
boundary between cases (C) and (A) occurs when the potential curves neither 
up nor down at the foot point s = (see figure i.e. when 

d^<i>,/ds%^o=0. (46) 

Assuming the foot point to rotate at the Keplerian rate fl = {GM/vjqY^^ , we 
find 

d^<i>e/ds^\s=o = n^isin^ 9-3 cos^ 9), (47) 

where 9 is the angle between the field line and the zu axis. Thus the boundary 
occurs at a critical angle (Blandford and Payne, 1982): 

9c = atan(\/3) = 60°. (48) 

The dependence on field line inclination is summarized in figure |[ On field 
lines more vertical than 60°, the situation is like in a stellar wind: there is a 
potential barrier to overcome, and this requires the existence of a hot atmosphere 
(temperatures comparable to the virial temperature). For inclinations less than 
60°, on the other hand, there is no impediment to the flow at all, and one 
would expect a large mass flux. This is, in fact, somewhat problematic, as 
discussed further in section ^ where the consequences of large mass fluxes are 
investigated. In between, there is a range in field inclinations (a narrow range 
if the disk atmosphere is cool compared with the virial temperature), where a 
reasonable mass flux results, and the magnetic wind theory works best. The 
detailed solutions of Konigl (1989) are of this type. It has been argued (Lubow et 
al. 1994) that the strong dependence of mass flux on inclination makes the flows 
in this range unstable. If this is the case, stationary magnetically accelerated 
flows may not exist, at least not from cool disks. The conditions for stationary 
flows may be better in AGN. Here, the likely presence of an ion supported torus 
(with ion temperature near the virial temperature, Rees et al. 1982) near the 
black hole would allow a magnetically generated flow from a wider range of field 
line inclinations. 

From this discussion, it will be clear that the details of how the wind is 
launched depend somewhat critically on things which are not presently known 
in detail, namely the inclination of the field lines near the disk surface and the 
presence or absence of a hot atmosphere. 

5.4 Geometry of the magnetic field near the disk 

In view of its importance for the wind launching problem, one would like to 
know what determines the shape of the field lines near the disk surface. Since 
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Figure 8: Launching conditions of a wind on field lines with varying inclination. 
The maximum of the effective potential along field lines is shown by a broken 
line. If the disk is cool, only a narrow range of inclinations around 60° yields 
'good' wind solutions. The slow, high rh flows at low inclination are problematic, 
see text. 



the field is close to a potential field near the surface, its strength and geometry 
in this region, including the inclination, is determined uniquely by a bound- 
ary condition at the disk surface^ namely the vertical component of the field 
strength. There are two possibilities for the origin of this field. If there is a 
dynamo process acting in the disk, one may expect field strengths of the order 
of equipartition with the gas pressure, as numerical simulations show (Hawley 
et al. 1995, Brandenburg et al. 1995). These simulations also indicate that the 
field is created with small length scales L, of the order of the disk thickness, 
in the radial direction, and somewhat longer in the azimuthal direction. The 
potential field created by such a small scale field decays with distance above 
the disk like exp(— z/L). This would not be the ideal field for driving magnetic 
winds. 

Another possibility is that the field is not internally generated, but is due to 
magnetic flux captured from the environment in during formation of the disk, 
and advected and compressed by the accretion process. Poloidal flux captured 
in this way cannot be destroyed by local processes in the disk, it can only escape 
by diffusing radially outward. The field strength would then be determined by 
the balance between outward diffusion and inward accretion (van Ballegooijen 

*A boundary condition at large distance is also needed. Since the field strength in the disk 
is likely to be much larger than in the interstellar medium, it is sufliicient to take a standard 
condition of vanishing field strength at infinity. 
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1989, Spruit 1994). In the absence of a theory for (turbulent) diffusion in an 
accretion disk, it is not possible to predict with any reliability what distribution 
of field strengths will result. Since all energy densities in the disk increase inward 
however, including that of the accretion flow, it is reasonable to assume that 
the balance will yield a field with inward increasing strength. The field above 
the disk will then have a shape like that suggested by figure ||. In Spruit et al. 
(1995) we have argued that fields of this configuration can be quite strong, with 
magnetic energy densities exceeding the gas pressure, which would make them 
ideal for the production of magnetic outflows. 



5.5 Poynting and kinetic energy fluxes 

The wind carries both kinetic and magnetic energy. The asymptotic ratio of 
these, at large distance, is a measure of how 'magnetic' the wind is. The Poynt- 
ing flux, in the MHD approximation, is 

S = — (vxB)xB. (49) 

The relevant component of S is that parallel to the poloidal flow. With ( p| ) 
and (p^), this can be written as 

S = -^n{w X B) X B • n, (50) 

An 

where n is a unit vector along Vp. Working out the cross products: 

s^nw^^. (51) 

47r 

Thus, the Poynting flux can be read as the work done by the rotation against the 
magnetic torque. At large distance, the azimuthal velocity (in the inertial frame) 
is small compared with ilw, so that by ( p^ Bp « B^Vp/^lw^ and 3> Bp. 
Thus 



where K is the kinetic energy flux ^pv^. Many flows (an example is the cold 
Weber Davis model, section |^) have their fast mode critical point at infinity so 
that (i'a/i'p)oo = 1, and q = 2. The magnetic and kinetic energy fluxes are then 
comparable at infinity. Near the disk, the Poynting flux dominates. Part of 
the Poynting flux is converted into a kinetic energy flux during the acceleration 
process. 

Relation (j5^) is valid only if the field survives in a highly wound-up form 
asymptotically. In section ^ I will argue that nonaxisymmetric instabilities 
are likely to destroy at least part of the toroidal field. In reality, the magnetic 
contribution to the energy flux may therefore be rather unimportant, and q <ti 1 
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rather than of order unity. If this is the case, we have the aestheticahy pleasing 
result that the magnetic acceleration process, after all its internal workings, 
produces a basically ballistic wind, which is only moderately magnetic. 

5.5.1 electron-positron flows 

There is some observational evidence for outflows containing electron-positron 
pairs (e*) from relativistic objects. In the galactic center source 1E1740.7- 
2942 a positron annihilation feature has been observed (Churazov et al. 1991, 
Churazov et al. 1994), and a transient feature has been seen in X-ray Nova 
Muscae (Gil'fanov et al. 1991). Thus there appear to be pair producing as 
well as jet producing black hole candidates (though no case is known yet which 
combines both), and it is natural to speculate that magnetic jets may exist that 
consist predominantly of a pair plasma. The Blandford-Znajek mechanism is 
thought to produce such flows (Blandford 1993 and references therein). Such a 
flow is technically not different from the jets considered above, since the same 
MHD equations apply. The main difference is that pairs may annihilate, thus 
removing mass and inertia from the flow. This would tend to increase the 
relative importance of the Poynting flux in the flow. At the same time, however, 
the Alfven speed would increase due to the decreasing mass density, and speed 
up the instability of the toroidal field. As the mass disappears, the magnetic 
field would therefore also disappear. In the nonrelativistic case, the flow would 
then decay completely into photons. If the flow is relativistic, it is possible that 
part of the energy of the decaying field is converted into a flux of low-frequency 
electromagnetic waves. For recent speculations on this topic, see Levinson and 
Blandford (1995). A magnetically driven pair plasma flow from a pulsar has 
been invoked by Arons and collaborators (Gallant and Arons 1994) for the 
Crab nebula. 

6 Cross-field balance 

The collimation of the flow is determined by the force balance in meridional 
planes, in the direction perpendicular to the field. The stationary, axisymmetric 
equation of motion that governs this balance is called the Grad-Shafranov or 
Grad-Schliitcr-Shafranov equation. For our case of a magnetized flow it has a 
somewhat complicated form. Some important aspects are discussed below, but 
for details I refer to Heinemann and Olbert (1978) and Sakurai (1985). To begin 
with, note that the solutions obtained in the above for a fixed poloidal field are 
still valid for the full problem, provided we read them as relations expressing 
the solution in terms of the (yet to be determined) poloidal field. 

The solutions of the azimuthal and longitudinal equations of motion are ( po[ ) 
[with (^] and (^. Inserting these into the original equation of motion (^, we 
get the required expression for the remaining, perpendicular, component. The 
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result is (Heincmann and Olbert 1978, Okamoto 1975, 1992): 
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(53) 



where a prime ' denotes d/dV'- It follows that the expression in braces must van- 
ish. This equation is to be read as a two-dimensional partial differential equation 
for the stream function z) of the poloidal field. Note that it is a somewhat 
implicit kind of equation, since it involves the quantities E{ip), K('ip),'cuAi'>P) 
which are known in terms of only as solutions of eqs. ( pO| ) and (^0|). Hidden 
in (^3|) is the fact that it is singular at the Alfven point. By working out the 
coefficient of the highest (second) derivatives, one finds that it vanishes at the 
Alfven point, so a regularity condition must be applied there. An additional 
complication in solving the equation is that it is of mixed type, namely elliptic 
in some parts of the (zu, z) space and hyperbolic in others. Where the bound- 
aries are is found out by computing the characteristics of equation This is 
conveniently done by inserting a short wave approximation: 



(54) 



and keeping only the highest (quadratic) terms in k. If /c|| and k±_ are the 
components of k parallel and perpendicular to Bp, this highest order treatment 
of the equation turns out to yield 
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K-«sp)(^'p -^^fp)' 
where a new critical velocity Vcp has appeared: 
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(55) 



(56) 



< V. 



As the fiow accelerates it first meets an elliptic region (fc^/fc^ < 0) for Vp 
then a hyperbolic region for Vcp < Vp < Ugp, another elliptic region v^p < Vp < 
Vfp, and finally another hyperbolic region for Vp > V{p. The significance of the 
critical velocity Vcp is seen by noting that it is the speed of the cusp of an 
axisymmetric slow mode wave (Heinemann and Olbert 1978). The surface on 
which Vp = Vcp is called the cusp surface. The cusp speed does not appear in 
the case of a prescribed poloidal field, since the wave mode involved bends the 
poloidal field lines. 

The characteristics of equation |5^ should not be confused with the character- 
istics of a time dependent MHD problem. Though various wave speeds appear, 
one is dealing with a time-independent flow. The ellipticity or hyperbolicity of 
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the problem refers to characteristics in the {w, z) space, not in an (r, t) space. 
PhysicaUy, however, there is a clear relation of the boundaries between elliptic 
and hyperbolic regions in the stationary problem on the one hand, and wave 
speeds in a time dependent problem on the other. This comes about because a 
wave in a frame comoving with the fluid appears as a stationary flow in the rest 
frame if the flow speed just cancels the propagation of the wave. This happens 
at the critical points. 

The presence of 4 different regions poses practical problems when construct- 
ing numerical solutions. The singular point at Vp = VAp has to be dealt with, 
as well as the boundaries between elliptic and hyperbolic regions, since different 
discretization schemes have to be used in each for numerical stability. The lower 
boundary condition at the disk surface, together with the regularity condition at 
the Alfven surface act as the boundary conditions for the elliptic regions (even 
though the boundaries of these regions do not coincide with these surfaces!). 
The solutions in the inner and outer hyperbolic regions are determined with 
values at the cusp and fast surfaces as initial data, respectively. If <^ Bp 
near the sonic point, the flrst hyperbolic region is quite narrow, and does not 
play an important role. One can then regard the entire region inside fast mag- 
netosonic surface as elliptic, also in practical solution algorithms. This is the 
case when the field is strong enough that 'interesting' degrees of acceleration 
take place. For discussions on numerical procedures, see Sakurai (1985, 1987) 
and Camenzind (1987). 



7 The character of the wind at high and low 

A simple model which demonstrates important parts of the physics is that of 
Weber and Davis (1967). I review it here in particular to discuss the dependence 
of the wind problem on the mass flux. This also relates to the question what 
happens to the flow when the inclination of the field line to the disk surface is 



less than 60° (cf. section 5.3 ) 



The model takes the poloidal field to be purely radial (in spherical coor- 
dinates), and looks only at the equatorial plane (with respect to the rotation 
axis)0. Though the model was invented for stellar winds, it can also be applied 
to the case of disk winds on low inclination field lines, nearly parallel to the disk 
surface. To further simplify the problem, I ignore the thermal pressure ('cold' 
limit), so that all acceleration is magnetic. 

Because he poloidal field is radial, its strength is given by 

Bp = Bo{wo/w)^, (57) 



®A purely radial field smaeks of monopoles. To remedy this, the field below the equator is 
given the opposite sign of the field above. The resulting configuration, with a current sheet at 
the equator, is physically realizable. For the dynamics of the wind, the sign of the magnetic 
field is unimportant. This is called the 'split monopole' configuration. 
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where wo is the foot point of the field fine on the disk, which we assume to 
rotate at the Keplerian rate O = {GM/w^Y^^. For the analysis it is practical 
to normalize nj and p to their values at the Alfvcn point, by introducing the 
variables 

x = m/zuA; y^p/pA, (58) 
and a normalized Bernoulli function 
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Substituting (p^) into (pO|), the Bernoulli equation then takes the form 
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Since the sound speed vanishes, tip = at the sonic point, and the gas density 
diverges there. With (|2^) the azimuthal velocity v'^ then vanishes also. Near the 
sonic point, the first two terms in (^^ describing the kinetic energy therefore 
vanish, and the only terms left are the gravitational and centrifugal ones. The 
condition dxH = then yields uj = xj^, in dimensional terms GM = fl^vu^. 
Thus, the sonic point Xs is at the foot point xo = wo/wa of the field line, and 



Xq 



UJ 



-1/3 



(62) 



The value of the Bernoulli function is now also known, 

1 1 n 3 



E 



^---ux, = --u 



1/3 



(63) 



The relation between (3 and uj follows from the fast point condition. It turns 
out that the fast point is at infinity (Goldreich and Julian 1970). We skip this 
part of the derivation. One finds then that x'^y remains finite at infinity, which 
corresponds to the fact that the fiow speed reaches a finite value at infinity. 
Using the conditions dyll = dxH = for the fast point, and the Bernoulli 
equation, all expanded for a; — > oo, the result is 



{x'y)f^'--uj-''\ 



/3 = -(i-c.-2/3)^ 



(64) 
(65) 
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Figure 9: The cold Weber-Davis model, showing Alfven radius wa = '^a/'^o, 
angular momentum flux J — J / [rj^^Q.ru'^) , asymptotic flow speed Voo/^'^Oi and 
field angle n — {B^/Bi-)a as functions of the mass flux. 



We need to express this result in more physical terms. For a given field 
strength and rotation rate, the external parameter determining the solution is 
the density at the base of the flow or, equivalently, the mass flux. Consider the 
mass flux as the external parameter. It is measured by the quantity r] = pVp/Bp 
[ci. (|l3|)], the mass flux 'per field line'. This has the dimension of the square 
root of a density. In fact, evaluating it at the Alfven point. 



The quantity 



Att) 
Bo 



(66) 



(67) 



Anflwo 

has the same dimension. It turns out to be the natural unit of mass flux in the 
model. It increases with the field strength, refiecting the fact that a stronger 
field is able to accelerate a larger mass flux to the same speed. Defining a 
dimcnsionless mass flux /i: 

Ai = 77/77,, (68) 

we want to express the results as functions of this dimcnsionless flux. From the 
definition of /?, and using (p2|), we find that 



Hence with (| 



(69) 



(70) 
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The location of the Alfvcn point is then 



mA/tno = c.i/3 = [^(l + M-'/3)]i/^ (71) 

When the mass flux is small, /i ^ 1, the Alfven radius is far from the origin of 
the flow. For large mass fluxes, /i ^ 1, toa does not get arbitrarily close to zuq, 
but reaches the minimum value 



t^o(3/2)i/2 (m^oo). (72) 



Further quantities of interest are, for example, the angular momentum flux per 
field line: 

j ^r]nwl = TjMujlfi^{l + fi-^^^). (73) 

This gives the angular momentum flux in terms of the mass flux and the condi- 
tions at the base of the flow. The terminal speed of the flow follows from ( p^ 
and (H): 

-^ = (/3c.)V6^^-V3. (74) 

This demonstrates one of the most important properties of the magnetic accel- 
eration model: it can produce wind speeds that exceed the escape speed flwo 
from the rotating object. In principle, it can accelerate a sufficiently small mass 
flux to arbitrarily large speeds, though in practice this ability is limited by the 



rather weak 1/3 power in (74) 



For 11 = 1, the final speed is just equal to the rotation velocity at the base 
of the wind, and for large mass flux, the final speed becomes arbitrarily small. 
What kind of flows are these massive but sluggish winds? A good way to see 
this is by evaluating the pitch angle of the field at the Alfven point. The model 
gives for this 

^EA = (2 -36^1/3 _/3 + ^)i/2^-i/2. (75) 

v^P/ A 



The limiting forms are 



p / A 



(19/8)1/2 (^«i) (76) 
1.14/i (/i>l). 



For small mass loss, the pitch angle of the field at the Alfven radius reaches 
a constant value which happens to be very nearly one radian. For large /Lt 
however, the azimuthal field dominates over the poloidal field. This is illustrated 
in figure |l^, showing the shape of the field lines for a large and a small value 
of /i. The case fj, ^ I can be properly called a centrifugally accelerated flow. 
Up to the Alfven radius, the field lines are not strongly bent, the fiow corotates 
approximately, and the poloidal flow speed can be found to good accuracy from 
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Figure 10: Shape of the field Hnes in a cold Weber-Davis model for low (left) 
and high (right) mass loss cases. Dashed line: Alfven radius. 

the effective potential. Beyond the Alfven radius corotation fails, so that the 
effective potential is not a good estimator for the flow speed any more. For high 
mass loss, the situation is very different. Corotation now fails right from the 
start, so that a strongly wound up (S^ » Bp) field develops long before the 
Alfven radius is reached. The flow is slowly 'pushed' outward by the pressure 
of the toroidal fleld, with final speeds much less than flrQ. Rather than being 
'flung out', the flow is more a sequence of magnetostatic equilibria, since the flow 
time scale w/vp is much longer than the dynamical time scales (GM/tn^)^^/^ 
and nj/wA- At the Alfven point {vp ~ WAp), for example, the ratio of the flow 
and Alfven time scales is vp^jvp = B^/Bp ~ fi. This disparity of time scales 
brings in the question of stability. If it is unstable, the highly wound-up field 
will change on the short Alfven time scale, interfering with its pushing activity. 
For the low-/i centrifugal case, the field is not strongly twisted, and stability 
is not an issue, until the flow reaches the Alfven radius. By then, most of the 
acceleration has already taken place. I return to the question of stability in 
section ^. 

7.1 Relativistic flows 

The relativistic case is somewhat outside the scope of this text. I will discuss a 
few basic properties of the special-relativistic case, and refer to Michel (1973), 
Goldreich and Julian (1970) and Li, Chiueh and Begelman (1992) for details, 
and Bekenstein and Oron (1978), Okamoto (1978,1992) and Camenzind (1987) 
for the general relativistic treatment. Also left out is the Blandford-Znajek 
model for magnetic flows driven by the rotation of a black hole. See Blandford 
(1993) and references therein. 

In the nonrelativistic case, the only parameter determining the behavior of 
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the flow was the dimensionless mass loss rate /x; the dependence on the other 
physical parameters could be found by simple scalings. In the special relativistic 
case, an additional parameter w = flzuo / c appears because the speed of light 
now fixes a velocity scale. Neglecting the gas pressure, a Bernoulli equation 
can again be derived as in section |^ (the derivation, in the this case, is most 
easily done in the inertial frame). It can, as before, be written in the form 
H(zu, p) = E, but analysing it is a bit more complicated. In the extreme- 
relativistic limit, in which the asymptotic Lorentz factor is large, the equivalent 
of relation (74) becomes 

loo = ^^-V3, (77) 
c 

(Michel, 1969), where fi is given by (^). As the mass flux is decreased, and 700 
increases, the Alfven radius asymptotically approaches the light cylinder radius 
c/r2. The high inertia of the flow at large Lorentz factor ensures that wa always 
stays smaller than c/fi. 

Eq. ( [ttI) shows that the flow can in principle become relativistic even when 
it is launched from a non-relativistic (flmo/c <C 1) object, if the mass flux is low 
enough. In practice, however, the weak dependence on /i means that Lorentz 
factors larger than a few can be produced only by relativistic objects. 

The cross fleld balance plays a more important role in the acceleration region 
than in the non-relativistic case. Whereas in the non-relativistic case the as- 
sumption of a flxed prescribed poloidal field is still fair near the Alfven surface, 
this is not the case for relativistic flows. In the extreme relativistic limit, the 
inertial forces in the flow are so high near the Alfven surface that they bend the 
poloidal fleld lines into a nearly horizontal shape at zua (Camenzind 1987). 



8 Collimation by 'hoop stress' 

The curvature force exerted by the toroidal fleld compresses the fleld conflgu- 
ration towards the axis. This effect becomes important only when the toroidal 
field is comparable to or larger than the poloidal field, otherwise the configu- 
ration is determined by the internal equilibrium of the poloidal field. For low 
mass loss flows (cf section collimation by the hoop stress in the toroidal fleld 
therefore starts roughly at the Alfven radius. The effect was flrst observed in 
calculations of the solar wind by Suess and Nerney (1975). 

Note that collimation here is meant in 'optical' sense: a flow is collimated if 
the flow lines are parallel. This says nothing about the width of the flow. For 
astrophysical jets, however, a collimated jet in practice is also narrow. This is 
because the central engine is very small compared with the scale of observed 
jets. An AGN jet for example that expands by a factor 1000 from its expected 
origin near the central black hole is still less than a parsec across. 

The asymptotic collimation of an initially radial flow is illustrated in flg- 
ure |ll|. In this model (Sakurai, 1985) a flow is launched spherically symmetric 
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Figure 11: Collimation of field lines in an initially spherical stellar wind. Radius 
scale is logarithmic, in units of the Alfven radius. All field lines asymptotically 
become parallel to the axis. The logarithmic distance scale distorts the field 
lines: on a linear scale the distance of each field line from the axis increases 
monotonically. From Sakurai (1985). 

on a 'split monopole' field. After passing through the Alfven point, the field 

becomes predominantly toroidal, and this causes the flow to become perfectly 
coUimated (all flow lines parallel to the axis) no matter how small the rotation 
rate of the star. The rate at which this collimation takes place, however, is slow, 
requiring several orders of magnitude in distance. For ordinary stellar rotation 
rates, the distance needed for full collimation is unrealistically large: the flow 
reaches its interstellar termination shock well before being collimated. 

In addition to the perfectly collimated flows like Sakurai's. solutions have 
been found that are asymptotically uncoUimated, in spite of the hoop stress. 
If the field configuration is such that the field lines diverge suSiciently rapidly 
near the Alfven surface (faster than a purely radial field) , the fast mode point is 
located at only a few Alfven radii, and the flow remains uncoUimated (Begelman 
and Li 1994). Depending on the conditions in the accelerating region, it seems 
one either gets an asymptotically fully collimated flow (the 'cylindrical' case) or 
an uncoUimated, space-filling flow (a 'conical' flow, see also Sauty and Tsinganos 
1994, Nitta 1994, Tomimatsu 1994, Heyvaerts and Norman 1996). In the latter 
case, the asymptotic ratio q of magnetic to kinetic energy flux is smaller than 
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in the cylindrical case [where q = 2, cf. section 5.5 1. 

The process of collimation by hoop stress is a natural consequence in ax- 
isymmetric rotating winds, and can be computed in detail. The completeness 
and accuracy suggested by such computations, however, is somewhat mislead- 
ing because they depend very heavily on the assumption of axisymmetry. If 
regions of predominantly toroidal field are as unstable as toroidal fields else- 
where in the universe and the laboratory, a significant revision of our picture of 
the collimation of magnetic winds is in order (section 5.5). 



9 Kink instability 

In the previous section we found that a predominantly toroidal field develops 
outside the Alfven surface. In high mass loss fiows, it develops also inside the 
Alfven radius. Consider first the case of a low-/i fiow, outside the Alfven surface. 
Assume that the flow is well collimated, and move into a frame comoving with 
the flow. In this frame, we see a toroidal field, slowly decreasing in time by 
the expansion of the flow. A predominantly toroidal field, however, is violently 
unstable to kink instabilities: such a configuration is equivalent to the linear 
pinch (e.g. Roberts 1967, Parker 1979, Bateman 1980). The mechanism of the 
instability is illustrated in figure ^ An initially axial, untwisted, magnetic field 
is wound up and becomes unstable when the azimuthal becomes larger than the 
axial field strength. This is akin to the instability of a twisted rubber band 
(figure p^). Instability sets in when the axial tension vanishes. Denoting by 
Bz and i?^ the axial and azimuthal components of the field, the axial component 
of the stress is {— B'^ + B"^) / Sir . The first term is the net magnetic tension due to 
the axial field, and is stabilizing; it likes to keep field lines straight. The second 
term, equal to the magnetic pressure exerted by the azimuthal component, is 
positive, expansive. When the pressure becomes larger than the tension, some 
of the energy put in by the twisting is released by a kink. Each kink reduces 
the number of windings by one, at the expense of increasing the energy in the 
axial field by lengthening axial field lines somewhat^. 

The kink instability is a transition to a nearby equilibrium of lower energy, 
i.e. the instability saturates at a finite amplitude. This is because the amount of 
azimuthal field energy that can be released is finite, while the energy expended 
lengthening the the axial field increases indefinitely with the amplitude of the 
perturbation. In a predominantly azimuthal field, the instability can also be 
visualized as shown in figure [T^. A stack of deformable disks (think of the 
disks in your spinal column, for example) is compressed (by the pressure of the 

^The condition > Bz can underestimate the degree of instability. A cylindrical field 
configuration typically becomes unstable already when it is twisted by more than one full 
turn, independent of the distance between the surfaces at which the twist is applied (Kruskal- 
Shafranov condition). In our case, this is not relevant, however, because B^/Bp increases 
with distance in such a way that the number of turns in the field is always less than one at 
the point where first exceeds Bp. 
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a b 

Figure 12: Sketch of the kink instabihty mechanism, a: in analogy with an 
overtwisted rubber band, b: when the azimuthal field dominates, the instability 
is like that of a stack of deformable disks under compression. 



azimuthal field). By slipping sideways and deforming somewhat, the disks can 
release some of the pressure, at the expense of the integrity of the stack. In 
stellar interiors, such configurations are also known to be highly unstable in 
spite of the presence of a stabilizing thermal buoyancy force (Tayler 1980, Pitts 
and Tayler 1985). 

In these descriptions of the instability, the field is treated as if embedded in 
a neutral medium, and the instability is a so-called 'external' kink. In practice, 
the flow could be surrounded by, or itself surround, a less twisted field configu- 
ration. This has a stabilizing effect. Conditions for instability in this case (the 
'internal' kink) are somewhat more complicated. For more on the subject see, 
e.g. Bateman (1980, Ch. 6)[[ 

Kink instability develops on a time scale ccj/ua^, the Alfven travel time across 
the flow, based on the azimuthal field strength. The moment that instability 
takes place, the azimuthal field providing the collimating hoop stress is reduced 
(see also Eichler, 1993). The energy involved goes into a less ordered field 
component which, if anything, adds outward magnetic pressure instead of an 
organized force towards the axis. Thus, the collimating effect of i?^ decays at 
the same rate as the instability takes place. The reduction of the collimating 
hoop stress has the strongest effect in the most collimated flows; this is seen as 
follows. If the coUimation angle is 6, the radial expansion speed of the jet is 
Vj^ = 9vp ~ Ova, which is small compared with the Alfven speed. Hence the 

'^In Bateman, and in the controlled fusion literature, the use of the words 'poloidal' and 
'toroidal' in case of cylindrical fields is opposite to our usage; this has to do with the torus 
geometry assumed there. 
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Figure 13: Production of a longitudinal field by kink instability. The longitu- 
dinal flow is taken to depend on distance from the jet axis. Initially parallel to 
magnetic surfaces, instability forces the flow to cross the displaced field lines. 
The differential flow speed stretches the displaced azimuthal field lines along 
the axis. 



instability has ample time to act as the jet moves outward. The effect of the 
instability would be less dramatic close to the Alfven radius. Choudhuri and 
Konigl (1986) have proposed that kink instability near the Alfven radius may 
be responsible for some of the alignment anomalies seen in jets at the VLBI 
scale. 

It takes longer than the instability time scale to dissipate the disorganized 
field component it produces (this is related to the known slow dissipation of 
magnetic helicity, and is seen also in numerical simulations, e.g. Galsgaard 1995). 
This dissipation, however, eventually leads to a reduction of the field strength 
compared with the standard axisymmetric jet. A second consequence of kink 
instability is therefore that the ratio of magnetic to kinetic energy flux in the 
jet becomes less than unity (see section 5.5). Since the Alfven speed is lower, 
the fast mode critical point is closer to the source, perhaps at only a few Alfven 
radii. Most of the observed jet would then be outside the fast mode point, 
and kinetic energy dominated. In short: the jet behaves like a ballistic flow, 
like a water jet from a fire hose. This would simplify the magnetic jet picture 
considerably: though the acceleration process is intensely magnetic, it would 
eventually produce a ballistically moving jet in which magnetic stresses are a 
secondary factor as far as the dynamics is concerned. 

Some observational evidence for the action of kink instabilities may be the 
fact that the magnetic field tends to be parallel to jet axis, at least in the 
faster (type II) jets (Bridle and Perley, 1984). If the flow speed along the 
jet is not exactly uniform over its cross section, the irregularities in the fleld 
produced by the instability will be stretched along the jet axis, see figure O. 
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The strength of this longitudinal field will be comparable to the kinetic energy 
of differential velocity across the jet. This field will have many small scale 
reversals of direction, explaining why the total poloidal magnetic flux inferred 
from observations (which are not sensitive to the direction of the field lines) is 
much larger than can be easily accomodated in the accelerating region. These 
observational indications can equally be explained by stretching of the field 
by interaction with an external medium, but irregularities produced internally 
by kinking have the advantage that they will also work in the absence of any 
interaction with the surroundings. 

10 Poloidal collimation 

In addition to 'toroidal' collimation by hoop stresses, a poloidal magnetic field 
surrounding the jet can also be a collimator. This is likely to be a powerful 
effect (Blandford 1993, Spruit 1994) if the disk is 'Large' (where by 'Large' I 
mean extending over a significant number of decades in radius). As an example 
to demonstrate this, assume that the vertical field strength at the surface of the 
disk is of the form 

^ {t'ItI + (78) 

where Ti is the inner edge of the disk and we take v to be between and 2. 
Then the field strength is largest in the inner parts of the disk, but the magnetic 
flux J rBzdr is dominated by the outer regions of the disk. This is a reasonable 
situation to expect for the field in a disk. A radially selfsimilar disk, for example 
(Blandford and Payne 1982), has v = 5/4. In the inner accelerating region of 
the flow such a field is close to the potential field given by the distribution 
of flux on the disk surface. It turns out that for a distribution like ([78[), the 
fleld lines have a nice, naturally coUimating, shape. An example is shown in 
flgure |lj (Spruit 1994), which shows the fleld lines for the case v = \. In this 
case, the field lines are parabolas, hence their collimation becomes perfect at 
large distance. Of course, the field stops to have this shape near the Alfven 
surface and beyond, and at distances where the finite size of the disk becomes 
noticeable. We can derive a maximum degree of collimation from these ideas, 
as follows (Spruit, Foglizzo and Stehle 1996). The best collimation is obtained 
when the Alfven surface is at a distance of the order of the disk size, but not 
further. At larger distances than this, the field more resembles that of a dipole, 
and does not have any coUimating properties. For a field of the form ( |78| ) one 
finds (for the case v = 1) that the angle of the field lines with the axis, at the 
Alfven surface, is 

0mi„ « {r-JraY'\ (79) 

where rj and are the inner and outer radii of the disk. If we assume that no 
further collimation takes place beyond the Alfven surface, for example because 
of the kink instabilities discussed above, this angle is also the minimum opening 
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Figure 14: Field lines of a potential field produced by a magnetic field strength 
varying as (r^ + on the disk surface. The coUimating shape of the field 

lines is due to the magnetic flux in the outer parts of the disk. 

angle of the jet. If the Alfven distance is significantly smaller or larger than rd, 
the collimation is worse. 

We can compare these minimum angles of collimation with conditions ex- 
pected for various kinds of disk. This is shown in table ||. This shows that 
poloidal collimation is capable of explaining opening angles of less than a de- 
gree in most systems, with the notable exception of Cataclysmic Variables. And, 
in fact, no CV is known to produce a jet, though there is evidence for outflows 
from many such objects (section I interpret this as a good case for the im- 
portance of poloidal collimation. A nice test case in this context is R Aqr. It 
consists of a white dwarf accreting from a giant companion, demonstrating that 
it is not the white dwarf nature of the primaries in CV that prevents them from 
having jets. Because of its very long orbital period, the disk in this system is 
probably several orders of magnitude larger than the disks in CV. 

Confirmation that the relative disk size rd/r; is important may perhaps be 
found in binary protostars. Several cases are now known (see Matthieu, 1996) 



Table 2: Typical disk dimensions and minimum collimation angle for poloidal 
collimation, for different kinds of accretion disk systems. 



protostars 0.01 AU 100 AU 

LMXB 10 km lO^km 

AGN 1 AU > lO'^Al 

CV lO-^km 2 lO^km 

R Aqr lO-^km > 10^ k 
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of relatively close binary protostars where at least one of the stars has a disk. 
The maximum size of such a disk can not be much larger than the tidal radius, 
something of the order of 1/3 of the orbital separation. The prediction is then 
that such disks do not produce well-coUimated jets if <, 30i?*, i.e. if the 
orbital separation is less than about 0.5AU. 
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